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Abstract 



Universal properties of the Coulomb interaction energy apply to all 
many-electron systems. Bounds on the exchange-correlation energy, in 
particular, are important for the construction of improved density func- 
tionals. Here we investigate one such universal property — the Lieb- 
Oxford lower bound — for ionic and molecular systems. In recent work 
[J. Chem. Phys. 127, 054106 (2007)], we observed that for atoms and 
electron liquids this bound may be substantially tightened. Calculations 
for a few ions and molecules suggested the same tendency, but were not 
conclusive due to the small number of systems considered. Here we extend 
that analysis to many different families of ions and molecules, and find 
that for these, too, the bound can be empirically tightened by a similar 
margin as for atoms and electron liquids. Tightening the Lieb-Oxford 
bound will have consequences for the performance of various approximate 
exchange-correlation functionals. 

1 Introduction 

The Lieb-Oxford lower bound [1] on the exchange-correlation (xc) energy is an 
exact quantum-mechanical property of any Coulomb-interacting system. Origi- 
nally, this bound was obtained in an investigation of a very fundamental problem 
of stability of matter [2] . An important practical application of the bound is as 
a constraint in the construction of approximate exchange-correlation functionals 
for density- functional theory (DFT). 

In DFT [Smi [S] , many-body effects on electronic structure are included by 
means of an approximate xc functional. The Lieb-Oxford bound is one of the 
key ingredients in the construction of nonempirical xc functionals [6J, where it 
is usually written as 




(1) 



where C is a positive universal constant whose exact value is unknown. A 
first estimate for the maximum value of this constant, Cl = 8.52 [2j was later 
substituted [T] by the tighter value Clo = 1-68, which is adopted in most current 
work. 

With this value of C, the bound is satisfied for instance by the local-density 
approximation (LDA) [3 [5], the generalized-gradient approximations (GGAs) 
PW91 [9] and PBE [10], and the TPSS meta-GGA [11]. On the other hand, 
earlier GGAs [T^] and meta-GGAs [13] and other functionals containing fitting 
parameters [TU [TSl [HI [H] are not guaranteed to satisfy the bound for all pos- 
sible densities. The success of some these latter functionals shows that, while 
the bound is doubtlessly obeyed by the exact functional, satisfaction is not a 
necessary condition for good performance in practice. 

In terms of the local-density approximation to the exchange energy. 



Clo = 1-68 is A^o = 2.275. The exact universal A and C are unknown, but 
their existence is guaranteed by the Lieb-Oxford bound, and their values are 
constrained by A < Alo and C < Clo, respectively. 

More recently, Chan and Handy numerically obtained an additional small 
reduction of this upper hmit of C to Cch = 1.6358, implying \ch — 2.215 
[T5] . The present authors [T^] empirically investigated the bound for atoms, 
model Hamiltonians and the electron liquid, as well as for a very restricted 
number of ions and molecules. It turned out that for all systems investigated 
there even the tightest previous formulation of the bound, employing Cch = 
1.6358, is still substantially too generous. Away from unphysical limiting cases 
of model Hamiltonians, Cconjec = 1-0 was found to be appropriate for all stable 
real systems investigated there. In view of the impact a reduction of C (or, 
equivalently, of A) may have on current density functionals, this observation 
requires further investigation. 

In the present paper, we enlarge the data base of the empirical investigation 
to include diverse families of ions and molecules, in order to accumulate further 
evidence and to be able to quantify the possible degree of tightening for molecu- 
lar systems. To this end, we use high-precision Quantum Monte Carlo (QMC), 
configuration interaction (CI) and exact exchange (EXX) data to evaluate the 
ratio 



for different families of ions and molecules, including first and second-row dimers, 
and various hydrides and oxides. 

This paper is organized as follows: in Sec. [2] we discuss various ions from the 
Be, Li, C and O isoelectronic series, as well as monovalent cations, comparing 






2 



1,22 
1,21 
1,20 
1.19 
1,18 
1,17 
1.16 
1,15 
1.14 
1,13 
1,12 
1,11 



O 



5K 

O X ■ 

o 



s 

o 

J6t O 



°oooooooooo 

o 



AAA 



A 



8 10 12 14 16 IS 20 22 

Atomic number Z 



Figure 1: Value of the ratio A, defined in Eq. (jj]), as a function of atomic number 
Z for different ions. Ions isoelectronic to He (diamonds), Be (plus), Li (circles), 
C (pentagons) and O (triangles), and monovalent cations (crosses) are shown. 
Values for neutral atoms (squares) |19j are included for comparison. 



their behaviour with data for the He isoelectronic series and atoms obtained 
in Ref. [19]. Section [3] deals with different classes of molecules, and Sec. |4] 
provides a comparative analysis of the emerging trends. Section [5] contains our 
conclusions. 



2 Ions 

Results from evaluating X[n] for ions from the He, Be, Li, C and O isoelectronic 
series, and for monovalent cations are displayed in Fig. [TJ 

According to Eq.Q, in order to calculate the parameter A[n] we need the ex- 
act exchange-correlation energy E^c [n] and the LDA exchange energy E^^^ [n] . 
Ionic correlation energies were taken from Ref. [2Dj, except for the He series, for 
which we employed the DMC xc energy given by Huang and Umrigar [3T]. To 
these near-exact correlation energies we added EXX energies calculated with the 
OEP method [551 US] to simulate the exact xc energy. Small differences between 
exact exchange in DFT and in Hartree-Fock theory are negligible on the scale 
of effects we are after. The values for E^^'^ were obtained 23| from the LDA 
exchange functional and evaluated at self-consistent LDA(VWN)[8] densities. 

Figured] shows that ionic A values are similar to those of atoms, and far away 
from the theoretical upper limits. X[n] decreases as the ionic charge increases 
and the density distribution becomes more compact. Chemically similar ions 
(same number of electrons) fall on common A(Z) curves, displaying systematic 
trends as a function of atomic number and total charge. Monovalent cations 
have charge distributions that are very similar to neutral atoms, in particular 
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for larger and consequently have almost the same A values. 

3 Molecules 

Figure reports A values for a set of dimers and hydrides. Figure [3] extends 
this comparison to oxides and a variety of other molecules for which near-exact 
correlation energies are available. 

For the first-row dimers we employed highly-accurate diffusion Monte Carlo 
(DMC) correlation energies from Ref. [25 {H2) and Ref. [25] (Li2 to F2). Cor- 
relation energies for the second-row dimers are from Ref. [26j . where they were 
extracted from experimental dissociation energies, corrected for zero-point os- 
cillations, by subtracting Hartree-Fock total energies. For the first-row hydrides 
we used DMC correlation energies extracted from Ref. [57], while correlation 
energies of second-row hydrides were obtained from experimental dissociation 
energies in Ref. [55]. For first and second-row oxides, as well as for the other 
molecules shown in Fig. [31 reference correlation energies are taken from [55] and 

El]- 

These correlation energies were combined with LDA exchange energies and 
Hartree-Fock exchange (again, small differences between exact exchange in DFT 
and in Hartree-Fock theory are negligible on the scale of effects we are after) 
to calculate the ratio A[n] = {E^'""=^ + Ef""""^)/ E^'^'^ . The Hartree-Fock and 
LDA calculations were performed with the Gaussian 03 package |29j. These 
calculations used the 6-311G-|-(3df,2p) basis sets for the first- row dimers and 
hydrides, and the cc-pVTZ basis sets for the second-row dimers and hydrides, 
and all oxides and remaining molecules. Geometries for the Hartree-Fock and 
LDA calculations were taken from the references providing the correlation en- 
ergies. For some molecules the geometries were not supplied together with the 
correlation energies [24j . In these cases we performed our own geometry opti- 
mization on the MP2 /cc-pVTZ level, followed by single-point calculations using 
the basis sets and methods indicated above. 

Note that in Figs. [2] and [3] we display results as a function of the number of 
electrons N of each molecule, instead of the atomic number Z used for atoms 
and ions. Nevertheless, the tendencies observed for ionic atoms clearly persist 
for neutral molecules: A(A'^) decreases for larger (heavier) systems and increases 
for more rarefied electronic distributions; chemically similar families of systems 
fall on common A(A^) curves; and values of \{N) are far from Xlo- 

The tendency of A to increase for more rarefied electronic distributions 
(lower average electron density) was already observed in the electron gas and for 
Hooke's atom [19]. To further investigate this tendency for molecules, we also 
considered the variation of A as a function of the distance between the nuclei 
of the H2 molecule. Approximate correlation energies of H2 as a function of 
bond length are available from Ref. [30] . From these data we find that for H2 at 
the equilibrium bond length \{Re = 1.401a.u.) — 1.230 (only slightly different 
from the value 1.251, obtained for the same system from DMC), whereas for 
stretched H2 this changes to \{Re = 5.0a.M.) = 1.409. Even this enhanced value 
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Figure 2: Energy ratio A as a function of the total number of electrons N, 
for first-row dimers (squares), second- row dimers (circles), first-row hydrides 
(triangles) and second-row hydrides (diamonds). 



of A, though, is still way below the maximum X^o = 2.275. 



4 Analysis 

The X[n] functional, introduced in Ref. [Hj, has two key properties. One is that 
by definition it measures the inverse weight of LDA exchange relative to exact 
exchange and correlation: The larger A[n], the smaller is |i?^''^^[n]| relative 
to |i?|^'^'^'[n]|. A[n] therefore serves to characterize the importance of nonlocal 
exchange and of correlation in a given system or class of systems, indicating 
when E^^-^l n] may be a good starting point and when more refined approaches 
are required. 

The other key property is that the maximum of X[n\ for a given class of 
systems provides a class-specific Lieb-Oxford bound. If these classes are large 
enough, such bounds may prove very useful in practice. Moreover, the maximum 
of A[n] across all investigated systems provides an empirical approximation to 
the universal constant A < Alo whose existence is guaranteed by the Lieb- 
Oxford bound ([3]), but whose exact value remains unknown. 

From the point of view of these two key properties, three common aspects 
of all A data presented in this paper deserve special attention: 

(i) A decreases as the density distribution becomes more compact {e.g., for 
increasing Z in isoelectronic series of ions). Larger A[n] arise for systems with 
spread-out rarefied density distribution. This trend is consistent with the char- 
acterization of very low-density systems as strongly correlated, which is rigor- 
ously correct for uniform electron liquids. 

(ii) Chemically similar systems are characterized by common trends of A, 
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Figure 3: Energy ratio A as a function of the total number of electrons N , for 
first and second-row oxides (crosses) and a set of miscellaneous other molecules 
for which near-exact E^c data are available. 



suggesting that A may be a useful parameter to characterize and classify sys- 
tems, and that A[n] may be a useful ingredient in the construction of refined xc 
functionals. 

(iii) All values of A obtained are much smaller than the maximum value 
XlOi averaging at around half of A^o — 2.275. The new ionic and molecular 
data reported here are fully consistent with the conclusion obtained in Ref. [19] , 
and reinforce the conjecture that the Lieb-Oxford bound can be substantially 
tightened. The empirical value suggested in Ref. [19 for real systems {i.e. 
excluding unrealistic extreme limits of model Hamiltonians), \conjec ~ 1-35 
(corresponding to C conjee ~ 1-0) still covers all new systems investigated here, 
for which A never exceeds 1.26. 

If we include unstable systems, such as H2 stretched to nearly four times its 
equilibrium distance, A[n] exceeds the value conjectured for stable systems at 
equilibrium (1.409 versus 1.35), which agrees with the trend observed for Hooke's 
atom in Ref. [19j. Even at equilibrium, H2 has with A = 1.251 the largest value 
found for any system investigated here, but this value is still clearly below the 
conjectured limit \eonjee ~ 1.35. 

5 Conclusions 

Numerical results for the Lieb-Oxford parameter A[n], based on high-quality 
input data for the correlation energy Ec [n] , confirm the expectation of Ref. [19j 
that for all realistic systems the upper limit Xlo — 2.275 is too generous, by 
almost a factor of two. This is the same conclusion obtained in Ref. [T^, but 
considerably strengthened by the present analysis covering many different types 
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of ions and molecules not considered there. 

Since we have not found a single exception to the emerging trend 31j, our 
analysis suggests to pursue the following further investigations: 

(i) Mathematical scrutiny of the Lieb-Oxford bound, which may not be the 
strongest possible bound for arbitrary Coulomb systems. Moreover, even if 
application to arbitrary systems {i.e. a truly universal bound) should require 
placing A at or near Alo = 2.275, "reaHstic" systems, away from unphysical 
limits of system parameters, may allow a much smaller value of the upper limit. 
It remains to be clarified how to characterize such realistic systems or parameter 
values in a way independent of A itself. 

(ii) Exploration of consequences a reduced value of A has for the perfor- 
mance of common approximate density functionals. This exploration is aided 
by an investigation of the emerging systematic behaviour of A[n] accross different 
families of molecules, to obtain analytical or numerical models of the functional 
AW. 

Work on these issues is in progress in our group. 
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